A STABILITY RESULT FOR THE STOKES-BOUSSINESQ 
EQUATIONS IN INFINITE 3D CHANNELS 
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Abstract. We consider the Stokes-Boussinesq (and the stationary Na- 
. vier-Stokes-Boussinesq) equations in a slanted, i.e. not aligned with the 

gravity's direction, 3d channel and with an arbitrary Rayleigh number. 
For the front-like initial data and under the no-slip boundary condition 
for the flow and no-flux boundary condition for the reactant temper- 
ature, we derive uniform estimates on the burning rate and the flow 
velocity, which can be interpreted as stability results for the laminar 
front. 



1. Introduction and the main results 

The Boussinesq system for a reactive flow is a model describing flame 
propagation in a gravitationally stratified medium [13]. It consists of the 
reaction-advection-diffusion equation for the reactant temperature (nor- 
malised so that < < 1) and the Navier-Stokes equations for the evolution 
of the incompressible flow u. The variables u and are coupled through the 
advection velocity in the reaction equation, and through the force term in 
the fluid equation. After passing to nondimensional variables [2], the sim- 
plified Stokes-Boussinesq system takes the form: 

(1.1) t + u-V0-A0 = f(0), 

(1.2) u t - vAu + Vp = Op, 

(1.3) divu = 0. 

Here, v > is the Prandtl number (inversely proportional to the Reynolds 
number), while the reaction rate is given by a nonnegative, nonlinear func- 
tion f(0) of ignition type. That is, / is Lipschitz continuous, and there is a 
threshold temperature < i?o < 1 such that: 
(1.4) 

f(0) = for < i? and > 1, f(0) > on (tf , 1), f'(l) < 0. 

The vector p = pg corresponds to the non-dimensional gravity g scaled 
by the Rayleigh number p, and we assume that p is non-parallel to the 
unbounded direction of the 3d channel D, which amounts to studying the 
system (1.1) - (1.3) in: 



D = (—00,00) xfi = {(x,x); x G R, x G fi}, 
1 
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with: 

P- e 3 + 0. 

The crossection Q is a sufficiently regular, connected and bounded domain 
in M 2 . For the solutions to (1.1) - (1-3), we impose the Neumann condition 
in 9, and the no-slip (Dirichlet) boundary condition in u: 

OA 

(1.5) — = and u = on dD. 

on 

The classical result by Kanel [7, 1] states that there exists the unique 
speed Co (necessarily positive) and a (necessarily decreasing) traveling wave 
profile <I> : R — > [0, 1], satisfying: 

-co$' - = /($), $(-oo) = l, $(+oo) = 0. 

We call Co and $ the laminar speed and laminar front, and to fix the ideas, 
we let 3>(0) = i?o- Then 9(x, t) = &(x — cot) is the unique (up to translations) 
traveling-wave solution to: 

(1.6) 9 t -e xx = f(6), 9(x,0) = 6 (x) 

joining the two equilibria 1 and 0. Moreover, if the initial data ^o^) differs 
from the Heaviside function H(x) by a compactly supported error, then 
there exists a shift xq such that, for the solution to (1.6) there holds: 

\\9{-,t) - $(• + c t + x )\\ L °°(R) ^0 as t ^ oo. 

Our purpose is to reproduce this result for the problem (1.1) - (1.3), (1.5). 
Following [3] and [2], we define the bulk burning rate B(t) and the Nusselt 
number N(t) by: 

B{ f) =ik I dx di ' B{t ) = 7 f B{ - s ) ds 

(1.7) ^ l Jd 1 h t 

N(t) = p J D \V6\ 2 dxdx, N(t) = j J N{s) ds, 

and the average flow U (t) by 

(1.8) U{t) = \ /Vc^loo ds. 

1 Jo 

Theorem 1.1. Assume that the initial temperature 0q(x,x) G [0,1] is such 
that 9q(x,x) — H(x) is compactly supported in D. Let uq G W 1,2 (D) and 
let (9,u,p) be a global solution of (1.1) - (1-3), (1.5) with: 9(-,0) = 9q and 
u(-,0) = uq. Then, there exists a constant Cq, depending only on Q,, such 
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o(Y)< B(l) < r„ f Co ( ^ + 4 ) 1 "(U- 



17 + Cn^r 
2 1^ 




(1.11) l/(t)<C n (^ + ^)+o(l). 



The above result shows that the solution of the initial-boundary problem 
for (1.1) - (1-3), with small p/u propagates with finite speed close to the 
laminar front speed. Also, note that if we replace 9 with the laminar front 
<1>, then a simple integration of (1.6) gives: B(t) = cq. This corresponds with 
the estimates in Theorem 1.1 when p = 0, i.e. the system (1.1) - (1-3) turns 
out to be a regular perturbation of the reaction-diffusion equation. 

Our next result states that with front-like initial data, the solution to the 
studied system stays front-like: 

Theorem 1.2. With the same assumptions as in Theorem 1.1, we have: 
(1.12) 

Cn,o 



$(x- cot + x + U(t)t + C a ,oVtj - < 9{x, x, t) 



Vt 

<<S>(x-c t-x + U(t)t + C^ov^) + ^ Vt » 1, 

with appropriate xo > and C^,o depending on £1, f and the initial data 9q. 
Finally, we have: 

Theorem 1.3. The results in Theorem 1.1 and Theorem 1.2 remain valid 
in either of the following cases: 

(i) the channel D = (—00, +00) x [0, A] is 2-dimensional 

(ii) the flow equation (1.2) is replaced by the stationary Navier- Stokes: 

—isAu + u ■ Vn + Vp = 9p, 

and the crossection £1 is sufficiently "thin", i.e it satisfies the same 
type of condition as in [10] : 

^Y^M 1/2 c P (c PW +(J\g. (o,x)| 2 d x y^ < 1. 

Existence of traveling fronts for the Navier-Stokes-Boussinesq system in 
infinite channels has been the subject of research during recent years [3, 2, 
11, 5, 8, 10]. Of interest has been also an understanding of the regulariz- 
ing and mixing effect of convection [6, 4]. In [3] and [2], the solutions of 
the system with front-like datum in a 2d strip have been considered and 
uniform estimates for the full Navier-Stokes-Boussinesq system have been 
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obtained for the stress-free boundary conditions on u. In this paper we gen- 
eralize these results to dimension 3 and the more physically relevant no-slip 
boundary conditions. The analysis follows [2] closely; it first seeks bounds 
for N and B using the parabolic equation (1.1). Where our argument di- 
verges from that of [2] is in finding an estimate for ||«||oo- I n [2], this has 
been done using Poincare's inequality for vorticity, based on the assumption 
that the vorticity vanishes at the boundary, and using the rather simple 
form of the vorticity equation in dimension two. In the present case, we rely 
on the almost-uniform Xie bound in Lemma 2.2. 



Acknowledgments. M.L. was partially supported by the NSF grants 
DMS-0707275 and DMS-0846996, and by the Polish MN grant N N201 
547438. M.R. was partially supported by the IPM grant no 90350025. 



2. Auxiliary results 

Our first result is a technical lemma, which extends Lemma 4.3 in [2], 
proved there for 2d channels (see also [4]). Here we prove a similar result in 
dimension 3, showing that diffusion in a tube behaves like Id heat equation. 

Lemma 2.1. Let u G W 1 ' 2 {D,W i ) with div u = be a given solenoidal 
flow satisfying u = on dD. Let <fi be the solution to the advection- diffusion 
equation: 

4>t + u ■ V(f) - A(f> = for (x, x) G D and t > 

(2.1) ^(x,x,t) = for(x,x)edD 

on 

4>{x,x,0) = (f> (x,x). 

Then there exists a constant Cq depending only on £1 (in particular inde- 
pendent of u and (J)q ) such that, for all sufficiently large t, there holds: 

II<K-,*)Hoo<^||0oHli(d) vt»i. 

Proof. 1. We first prove a Nash-type inequality, valid for solutions of (2.1): 



(2.2) iiv<x,t)in 2(D) > c a j^rrm iu? ,,113 v *- 



To simplify the notation, in what follows we suppress the dependence on t 
and write <p instead of <f>(-,t), etc. 

Define ip{ x ) = J^ < X a; i^) dx and let k(x,x) = (f>(x,x) — ip{x), so that 
4> = ip + k and j n k(x, x) dx = for every Notice that, by Cauchy- 

Schwarz inequality: 

(2.3) |M|| 2(D) = \n\ jf (£ df) 2 dx < jf jT |^| 2 = u\\l HD) , 
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Similarly: 



, . 11^ II 1,2 (D) < \\<Px\\l2(D), IIV'llLi(D) < \\P\\Li(D), 

\\k\\ L 2 < 2\\cf>\\ L2 , \\Vk\\ L 2 < 2||V0|| L2 , ||fc|| L i < 2\\<(>\\ L i. 

Let V> : R — > C be the Fourier transform of ?/>, i.e. ^(w) = ip{s)e~ 2muls ds. 
By Plancherel's identity, we have: 

( 2 - 5 ) II^IIl2(R) = IHIl2(R); ll^'llz^ffi) = ||27TK^||| 2(ffi) . 

For a given positive m, we now write: 



./|w|<m ^|w|>m 

<2mWi^) + i / w 2 |^( W )| 2 d^ 
<2m||^||i 1(R) + i ^||V , |li 2(R) 

where the estimate of the first term follows by: |Y>(w)| < H^Hl^d); while to 
estimate the second term we used (2.5). 
Setting m = II^II^rjIMIIi 2 ^). we obtain: 

wi> w <(^ + 2)11^11^11^11^, 

which implies: 

Mho) < (^+2)i^r 2/3 iiV'iil( 3 (D) ii^ii 2 : / 2 3 (D) - 

Now, by Cauchy-Schwarz inequality, the Sobolev embedding W 1,2 (D) ■— > 
L 4 (D), and the Poincare-Wirtinger inequality on if follows that: 

\\k\\l HD) < \\k\\f {D) \\k\\% 3 {D) < Cn\\k\\f {D) \\Vk\\% 3 {D) . 
Therefore, by (2.3) and (2.4): 

Mh (D) <2ml HD) +nk\\h (D) ) 

< cn (u\\f {D) \\m% 3 {D} + \m¥ {D) \m\%] D) ) ■ 

We now argue as in [2]. Since y := ||V0||^2 3 satisfies: ay 2 + by — c > 
with appropriate a, b, c > 0, then: y > =^±^ = —0== > 
Hence: 

iiv^ 3 > ihi^ 

which gives: 
yielding exactly (2.2). 
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2. Recall now that: 

( 2 - 6 ) IHIli(£>) ^ II^oIIl^d)- 

Indeed, the L 1 norm of a solution to (2.1) is conserved when the initial 
data is positive. In the general case one can write (po = (pQ — (pQ where 
4>q and 4>q are positive, with disjoint supports. Solving (2.1) for each, one 
obtains the inequality (2.6). Further, integrating (2.1) against (j) and using 
incompressibility of u and the boundary condition, it follows that: 

(2.7) ± m i 2 = -2\\vm 2 . 

In view of (2.2), (2.6), (2.7) we now obtain: 



d* ^~"\\<h\\ii + \\<hh4<i>\\h 
which, after integrating in time, gives: 

(2.8) t<Cn(^ + U ° hl 



L 2 \m\L* 

Call a = IMIWH^oIIli- From ^ 2 - 8 ) il follows that ^TT - The function 

4 

a I—?- is increasing and it converges to as a — > 0. Let j3 be the unique 
solution to J^ +1 = ^p, so that a < (3. Now, for t — > oo clearly ^jp- — >■ 0, 



hence also (3 — > and 4>- < = ^p. Consequently, f3 < and we 



arrive at: 

Ml* < %\\MW Vt»l. 

3. We now argue as in [2]. Let Vt be the solution operator for (2.1). Thus 
far, we have showed that: 

\\Vt\\ L ^v<^ Vt»l. 

Now if VI is the adjoint operator, then V\ is the solution operator to (2.1), 
with u replaced by —u. Therefore the above argument works again: 

\\n\\ L ^<% vt»i. 

Finally, we conclude the lemma: 

= \\Vt\\L^LA\V*t\\L^ < || vt » i. 



We now present a lemma taken from [10]: 
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Lemma 2.2. Let g G L 2 (D). There exists a constant Cq,, depending only 
on the crossection VL, such that for any solenoidal flow u G W 2 ' 2 n W ' (D) 
satisfying: 

—vAu + Vp = g, div u = in D, 

there holds: 



(2.9) IMIoc < ■^=\\Vu\\^ 2 2 {D) \\g\\ 1 / 2 2 (D) + Cn||Vu|| L2(D) . 

We remark that the proof of (2.9) relies on Xie's estimate [12]: 
\\u\\ < _L||Adl 1/2 HVdl 1/2 

\/27T L2(D) L2(D) 

valid for u G W 2,2 (D) n Wq' 2 (D), and on a commutator estimate [9]: 
|| (PA - AV)u\\ 2 LHD) < Q + ||Au||| 2(D) + C Ae ||Vu||| a(JJ) 

where V is the Helmholz projection onto the space of solenoidal vector fields. 
The interest in the inequality (2.9) lies in the independence of the constant 
at the term involving g. Indeed, this was the key argument allowing to 

prove [10] existence of traveling waves for the full Navier-Stokes-Boussinesq 
system in 3d channels, satisfying appropriate "thinness" condition on the 
crossection £1. The same argument is needed to obtain the uniform bounds 
for the stationary Navier-Stokes-Boussinesq system in Theorem 1.3 (ii). 

On the other hand, by elliptic estimates and the Sobolev imbedding, it 
follows directly that: 

(2.10) IHloo < Cn\\u\\ W 2,2( D) < Cn (\\^u\\l 2 {d) + ^Wah'^D)^ ■ 

In fact, already this inequality is sufficient for the estimates in case of the 
Stokes-Boussinesq system. 

We will also need the following result, in the line of Lemma 3.6 from [10]: 
Lemma 2.3. For each t, there exists a function h 6 W^ 2 (D), such that: 

\\e(;t)p-vh\\ L 2 {D) < c n p\\ve(;t)\\ L 2 {D) , 

with Cq depending only on Q. In fact: 

(2.11) Cn = C PW + (£|5-(0,x)| 2 dx^j 1 , 

where Cpw stands for the Poincare-Wirtinger constant of Q. 

Proof. Let ei,e2,es be the standard basis for M 3 . Suppressing the time 
variable t, we define: 

h(x, x) = p ■ e\ \ -r 9(x, x) dx dx + p ■ (0, x) -r 9(x, x) dx. 
Jo Jn Jn 
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Therefore, the following identity concludes the proof of lemma: 

Op — V/i = p I 9(x, x) — -f 0(x, x) dx I — p • (0, x) -r 9 x {x, x) dx e\, 
\ Jn J Jn 

by the Poincare-Wirtinger inequality: \\9 — j 9\\L?(n) < Cpw||V0||^2^. ■ 

3. Proofs of the main result 

The following lemma has been proven in [2] for the case of 2d channels 
and vorticity of the flow u vanishing at dD. Exactly the same proof, relying 
on the construction of super and sub-solutions to (1.1) is valid also in the 
present 3d case. Since the argument uses the estimate in Lemma 2.1, we 
partially reproduce it below for the sake of completeness. 

Lemma 3.1. There exists a constant Cap depending on 17, /, and on the 
initial data 6q, such that: 

(3.1) N(t) < l -B(t) + U(t) + C Qfi (j + V* » 1. 

Moreover, there exists xq > and q G L 1 (]R) ; such that for all sufficiently 
large t » 1: 

$ (x - cot + x + U(t)t + Cn j0 Vt) - Q(x, x, t) < 9{x, x, t) 
(3-2) V 1 _ 

< $ ( x - cot - x - U(t)t - CnfiVtj + Q(x, x, t), 

where Q is the solution to: 

Qt + u ■ VQ - AQ = in D, 

(3-3) 3Q 

-— = on dD, Q(x,x,0) = q(x). 
on 

Proof. We only prove (3.2), since it implies (3.1) as in [2], Lemma 4.2. 
Define: 

tpi(x, x, t) = $ (x - cot + x + U (t)t + C\/t) - Q(x, x, t), 

where xq,C > are to be determined later, and Q is as in (3.3) with q 
appropriately chosen. 

To prove that tp is a subsolution, we first need to show the non-positivity 
of the following expression: 

(3.4) 

{^i) t + u-V^i- Atpi - f(ifi) 

t)h°°(D) + + ni(x,x,t)^j $' (x - cot + xo + U(t)t + CVt) 

+ /($)-/($-Q) 



STOKES-BOUSSINESQ EQUATION IN A CHANNEL 9 

Take q G L l (R) such that < q(x) < a = min(tf /2, (1 - #o)/4). By the 
maximum principle we have: 

(3.5) < Q(x,x,t) < a. 
Let now xq > be such that: 

(3.6) 9 (x,x,0) > $(x + x ) - q(x). 

Finally, let C be large enough so that when $ G (i?o, 1 — (1 — #o)/4) then: 

^=d>' + /(*) - /(* - Q) < ^=d>' + ||/'|| L oc \\Q(; t)|| £ « (D) 

where we used Lemma 2.1 to estimate ||Q||oo- On the other hand, for $ G 
(0, i?o) U (1 — (1 — #o)/4, 1) the nonpositivity of the right hand side in (3.4) 
follows directly, via (3.5). Concluding, (3.4) and (3.6) imply the lower bound 
in (3.2). 

The upper bound follows by assuring that: 

1p r (x, X,t) = $ (x — Cot — Xq — U(t)t - Cy/tj + Q(x, X, t) 

is a supersolution. Similarly as above, this follows by choosing xo such that: 

6q(x, x, 0) < $(x - x ) + q{x) 
in addition to (3.6), and having C large enough so that: 

(W)t + u-VW - AW - /(W) > -^jf' - ||/'||l«»||Q(-,*)|| l =o (d) 

>-^7= min |$'( s )| -^M||/'|| Loo ||g|| > , 

when$€(0 o /2,l-(l-0o)/2). ■ 
Clearly, Theorem 1.2 follows from (3.2) and Lemma 2.1. We now have: 

Lemma 3.2. There exists a constant Cq^ depending on Q, f, and on the 
initial data 6$, such that for all sufficiently large t » 1: 

(3.7) B(t) < c + U(t) + C n ,o Q + ^) > 
and 

(3.8) B(t) > co - tf(t) - C7n,o Q + ^) ■ 
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Proof. Denoting = U (t)t + Cn ) ov / iand using the bound (3.2), we obtain: 
(3.9) 

= / 9(x, x, t) — 9(x, x, 0) dx dx 
|"F Jd 

< j^t- I $(x - cot -x - £(t)) + Q(x, x, t) - $(x + x ) + Q(x, x, 0) dx dx. 
|"|* Jd 

By (2.6) and the construction of the corrector Q, it follows that: 

(3.10) 77777/ \Q(x,x,t) + Q(x,x,0)\ dx dx < -\\q\\ L itm\ < 
\il\t Jd t 

Further: 



- - x o ~ £(*)) ~ + x o) 
= T / - c t - so - £(*)) - 1) 

t J -oo 

/-xo /•cot+xo+C(t) 
(1 - $(x + X )) + / (1 - $(x + X )) 

-oo J —xo 

roc 

(3.11) + / $(x - c t - x - £(t)) 

•/c t+x +£(t) 
/ — xo /-OO 

- / 3>(x + x ) - / $(x + x c 

J c t+x () +£(t) J —xo 

i r^o+m ^ ^ | Co t + 2x + g(t)| 
< c + £7(t) + C n ,o (j + , 

which together with (3.10) proves (3.7). To prove (3.8), similarly as in (3.9) 
we note that: 

B(t) > — — / $(x-cot+x +^(t))-Q(x,x,t)-$(x-x )-Q(x, x,0) dx dx, 

and as in (3.11) we obtain: 

1 , , , ™ n |c t-2x -^(t)| 



/ $(x - c t + x + £(t)) - $(x - x ) 
> c - U{t) - C n , ( 7 + 



\n\t 



t Vt 

Together with (3.10) the above implies (3.8) 
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Proof of Theorem 1.1. 

Multiplying the fluid equation (1.2) by u and integrating over D, one obtains: 

d o 2 

(3.12) -dMIia + v\\ Vu \\h < Cn—\\V0\\h- 

at v 

Integrating (1.2) against ut gives, in turn: 

\\ u t\\h + < Cnp 2 \\ve\\ 2 L2 , 

where in both inequalities above we used Lemma 2.3 to "replace" the term 
Op by pVO, Taking averages in time, we get: 

ft J2 



(3.13) 



7 f ||Vu||£ 2 dt < C n ^N(t) + -\\Vu \\l„ 

t Jq V Ut 

7 f \\u t \\ 2 L 2 dt < C nP 2 N(t) + 7||V« ||1 2 - 
1 Jo 1 



By Lemma 2.2 or (2.10), and Lemma 2.3 it now follows that: 
IMIoc < C n M|Vu|| L 2 + -\\u t \\ L 2 + ^||W|| L2 J . 



Using (3.13) we obtain: 
(3.14) 

m < ((i jT iivniii.) 1/2 + 1 (i £ i K n| 2 ) 1/2 + ^ 

< c7 n (^y/Wyj + -^=||Vu || L 2 + -L||v«olUa) 

By (3.1) and (3.7) we hence get, for large i >> 1: 

*w<i* + §ew + Gv,(i + ^) 

(3.15) 

where C a // depends on £l,f,v,9o and uq. Consequently: 



and, returning to (3.14): 



In view of Lemma 3.2, Theorem 1.1 is hence proven. 
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Proof of Theorem 1.3. 

We only prove the assertion (ii), because the 2d case in (i) follows with 
exactly the same calculations as in Theorem 1.1 and Theorem 1.2 . 

For the stationary Navier-Stokes-Boussinesq system, using the Poincare 
inequality and (2.11), we obtain the following counterpart of (3.12): 



where Cp and Cpw are, respectively, the Poincare and the Poincare- Wirtinger 
constants of Q. By Lemma 2.1 and argueing as (3.14), we arrive at: 



It is therefore clear that when the constant in front of N(t) in the last term 
of the right hand side above is smaller than 1, the results of Theorem 1.1 
and Theorem 1.2 follow as in the case of the Stokes-Boussinesq system. ■ 
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